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Abstract. We investigate the influence of an interband exchange interaction on 
magnetism in a two-band Hubbard model. Our main emphasis lies on spin-flip 
scattering which is often neglected but is neccessary to retain the full rotational 
symmetry of the Hamiltonian. We find a striking dependence of the magnetization 
on the interband exchange coupling constant J and a substantial suppression of 
ferromagnetic order for a large range of values of J . The onset of an RKKY-like 
magnetic ordering mechanism is also observed. 



PACS numbers: 71.10.Fd, 71.20.Be, 75.10.Lp 
1. Introduction 

Electronic correlations in the transition metals Fe, Co, Ni still present a major challenge 
in condensed matter physics. The stability of ferromagnetism in these materials is not 
yet fully understood. Contrary to rare-earth systems, these materials have no localized 
magnetic moments on "atomic" orbitals. The formation of finite magnetic moments in 
an itinerant electron system has to be explained. 

In the past, there have been many attempts to describe transition metals 
theoretically ||, H, H, ||, H, 0]. Virtually all of these are based on a modification 
of the multi-band Hubbard model @, |^ , in which only the on-site Coulomb interaction 
is considered. Recently, the importance of band degeneracy for correctly describing 
ferromagnetism has been confirmed by Quantum Monte Carlo (QMC) calculations [0, 



TT| . But since all calculations mentioned above rely on some approximations (even 
the quasi-exact QMC calculation works only on a simplified Hamiltonian), no complete 
understanding of the complicated d-band metals is reached yet. The only way to gain 
a picture of the complex physics of these systems is to consider partial problems for 
which conclusions can be drawn from the available approximate approaches to more or 
less oversimplified theoretical models. In this paper, we want to provide another small 
piece belonging to the puzzle of the (i-band metals. 

The investigation of ferromagnetism in the Hubbard model has a long history. 
In fact, it was the original intention introducing this model It is known 

for a long time that for a band filling of one electron above (above or below) half 
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filling in the limit U ^ oo, the ground state is ferromagnetic for a fee- (sc- or bcc- 
) lattice(Nagaoka state) p3| , The stability of the Nagaoka state was subject to 
extensive investigations. E.g. in infinite dimensions {d = oo), its stability could 
be proven for a wider range of band fillings and finite U |]T5[. Using variational 
treatments, further limits could be set on its stability on various lattice structures in 
two or three dimensions [T^. Related to these statements is the so-called fiat-band 



magnetism ITT]. Here, a ferromagnetic ground state could be rigorously proven for 
a dispersionless band structure. Furthermore, for finite temperatures, the existence of a 
ferromagnetic phase for certain parameter ranges has been established using dynamical 
mean-field theory |19] . We therefore believe that the intraband Coulomb interaction 
as described by the single-band Hubbard model is one major ingredient for itinerant 
ferromagnetism. 

However, besides the strong on-site Coulomb interaction, transition metal systems 
are also characterized by the 5-fold degeneracy of the d-bands. How does this fact 
influence ferromagnetism? In atoms, Hund's rules will favor a parallel alignment of 
the spins of electrons on degenerate levels. That orbital degeneracy will enhance 
ferromagnetic stability also in lattice systems appears to be a fact. The general validity 
of Hund's rule for the ground state of a degenerate Hubbard model has indeed been 
proven (for further statements concerning ferromagnetism in degenerate Hubbard 
models see pl| , 0). Often, this is the justification for using a simplified interband 



exchange interaction, which is restricted to a longitudinal Ising-like spin exchange. 
The full SU(2)-symmetric interband exchange interaction can be separated into a 
longitudinal (Ising) and a transversal (spin-flip) term as will be shown below (see 
equation (|^)). The longitudinal interaction will try to align the z-components of the 
spins due to the energy gain connected with a positive value of interband exchange 
constant J, thus fulfilling the predictions made by the analogy to the atomic behaviour. 
This is most obvious in mean-field theory, where the transverse part vanishes, and a 
magnetization-dependent bandshift is induced by the longitudinal component. 

In several recently published papers, more sophisticated calculations for multi-band 



Hubbard models were presented, as e.g. in Quantum Monte Carlo [10, 22], slave- 



boson p3| or Gutzwiller- variational ansatz methods [^. But many of these neglected 
the transverse part of the spin exchange with the reasoning explained above. Up to now 
it is not clear what the influences of the disregarded terms are. 

A valuable contribution to the problem of ferromagnetism in orbitally degenerate 
Hubbard models was given in where the authors use an exact diagonalization 
method on a restricted Hilbert space of an in the limit of d = oo equivalent impurity 
model. 

The aim of this paper is to investigate the influence of spin-exchange processes 
originating from the transverse part of the interband interaction on ferromagnetic 
stability. We examine a minimal two-band model which includes those parts of the 
general Coulomb interaction that we believe to be the most important for the stability 
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of ferromagnetism [Q, ^, ^, 0]. This leads to the following Hamiltonian: 

H = Ho + Hu + Hj (1) 

Ho= ^™(^) "fc,m,aafc,m,a (2) 

k,m,cr 

Hu = ^ ni^rn,cT ni^rn-<T (3) 



"^•^X^ '''j.Tt • CTj.m (4) 



We use the usual notation for the electron annihilation (creation) operators CL^^^n a 
with wavevector fc, band index m (m being the complementary band of the two-band 
system) and spin o. The free bands are described by the dispersion em(fc), U and J 
are combinations of the appropriate Coulomb matrix elements as introduced e.g. in 
0, H, 1^ . The spin operators in Hj are defined as 

1 V- 



with z„ = +1 (—1) for spin | (spin j). Within this Hubbard- type Hamiltonian, the 



intraband part Hu is able to produce ferromagnetism for sufficiently large U 0, [26 
The second interaction, Hj introduces interband exchange processes of two different 
kinds, as already mentioned above. This can be seen in the following decomposition of 
Hj 

Hj = --J ^ {cFlm(^i,^ + Cflm(^i,m) (6) 

The first term (transverse part) represents spin-fiip scattering and the second 
one (longitudinal) an "Ising-like" exchange, which tends to stabilize spontaneous 
ferromagnetic order. Many calculations on multi-band Hubbard models only consider 



the Ising term and neglect the spin-fiip part of the interaction [|T0], ^ |23|, |2^. In this 
paper, we want to trace both parts of the interband interaction (^ with identical quality. 

It goes without saying that (H) and (|) do not respresent the full set of local 
Coulomb interactions between (i-electrons. Our previous studies 0, ^, 0, however, 
have evidenced that they are most important for treating magnetic phenomena in 
transition metals. The final goal of our investigation will be the more or less quantitative 



description of real substances such as Fe 0, Co |]^, Ni g], and Gd ||2S[|. For this purpose 
we combine many-body model methods with "a6-^n^^^o" -bandstructure calculations. By 
definition the underlying many-body model incorporates only those interactions which 
are believed to be decisive for the collective magnetism with respect to temperature 
dependencies and typical correlation effects, and which are probably not properly taken 
into account by usual LSDA-treatments. According to our previous calculations the 
interactions (|) and (^ should be most important while, e.g., the interband Coulomb 
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interaction (denoted by U in our previous studies 0) turns out to be not that decisive 
for the band-ferromagnetism of transition metals. It can be assumed that this part 
of the local Coulomb interaction is sufficiently well covered by an LSDA calculation. 
In this sense our present model investigation is to be understood. Instead of tackling 
a more complicated, but probably still insufficient Hamiltonian, we aim at a better 
understanding of those interactions which might be the roots of the phenomenon 
ferromagnetism, without referring to a real system. The description of more realistic 
systems seems possible by a combination of our model study with LDA calculations 
along the kines of ref. p|, ||, 0], which is intended for the future. 

In the following section, we introduce an effective medium method to decompose 
the complicated many-body problem of the Hamiltonian (|I]) into two separately solvable 
problems of simpler structure. After that we introduce the approximations that lead us 
to a fully self-consistent solution for the two-band model. In section ||, we present and 
discuss our results. 

2. Theory 

2.1. Effective Medium Approach 

Even though we restricted ourselves to a rather simple model Hamiltonian, one is in 
need of a convincing approximation method to solve the problem. We propose a method 
based on an effective medium ansatz. This ansatz will map the original problem of 
Hamiltonian (P onto a set of simpler model Hamiltonians, for which well tested standard 
approximations exist. These can be put together in an appropriate way to get a solution 
of the original problem. This method is generalizable to many different models, a similar 
approach applied to the periodic Anderson model was recently published [p^ ]. 

We want to introduce the method using the two-band Hamiltonian ([^) with the 
two interaction terms Hu and Hj. The self-energy Tjk,m,cr{E) can be defined using the 
equation of motion of the single-electron Green's function Gk,m,a{.E) = {{ak^m,a] (^kma))'- 

EGk,mAE) = fi + em{k)Gk,mAE) + (([ow, Hu + Hj]_-al,^J 

= h + tm{k)Gk,m,a{^) + '^k,m,aiE)G k,m,aiE) (7) 

Using the linearity of the commutator and the Green's function, one can define "self- 
energy parts" , 

((K^,., HuU al^J = T.%lAE)Gu,mAE) (8) 

((K^,., HjU al^J = ^ulAE)Gk,mAE) (9) 

with T^h^mAE) = Sfc2i cr(-^) + ^fc'm ui^)- Assuming the knowledge of either one of these 
self-energy parts, one can introduce the following effective Hamiltonians: 

^'"^ = E (^-(^) + ^uIM) + Hu (10) 

k 

H^!r^ = Yl i'rnik) + J:^^iAv)) <m,.ak,m,. + Hj (11) 
k 
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which formally depend on a parameter r]. By solving each of these Hamiltonians for all 
values of rj, one can obtain the respectively missing self-energy part using the following 
identities: 

ss,.(i?)=4:s^Hi?)i.^i. (12) 

^ZAE) = Ei^ym,=E (13) 

with ^k^f\'E) bsing the self-energy of the effective Hamiltonian ([TUp and '^k^'^a\-E) 
of Hamiltonian (|11]). The introduction of the energy-parameter i] in the effective 
Hamiltonians (|TUp and (^) is neccessary in order to distinguish two different kinds 



of energy. The effective medium energy parameter r] should not be confused with the 
energy as dynamic variable of e.g. an equation of motion method to solve the respective 
partial many-body problem. This distinction is mandatory since otherwise, the system 
could show unphysical behaviour. Furthermore, for the same reasoning it is neccessary 
to calculate expectation values in the full system, e.g. by using the formal solution of 
equation (J^) for the single-electron Green's function: 

G^,„.AE) = 7 7^ rj, ^ (14) 

The two problems posed by the Hamiltonians ( ]T0| ) and (|1T]) are strongly related to each 
other. The solution of one in form of the respective self-energy (^) or (|T^) is needed 
as input for the other. This implies a self-consistency condition on the two self-energy 
parts which can only be fullfilled in an iterative way. But the advantage of the effective 
medium ansatz is also rather obvious: The two Hamiltonians (|10D and (|ll]) have been 
analyzed very well by now since they are both standard models of many-body-theory. 
Model (|l^) is essentially the single-band Hubbard model , and model (|ll|) is known as 



s/-model or ferromagnetic Kondo-lattice I^O], ^ . So for both partial problems already 
known approximation schemes can be used. 

For the Hubbard-model, many useful approximations exist. We use the spectral 



density approximation (SDA) ||26|, ^ which has to be considered as a strong- 
coupling theory. Although this choice prohibits an investigation of the small-t/ 
behaviour, its advantages are enormous. Besides of being mathematical simple and 
numerically reasonably fast, it has proven to give a qualitatively correct picture of 
ferromagnetism in the strong-coupling regime and compares to more sophisticated 



approaches It will be described in the following chapter. The s/-model has also 
attracted much interest. A very promising interpolating, moment-conserving equation- 
of-motion decoupling scheme has been developed in refs. |^ ^ . This method will 



be discussed in section p.3| . The combination of these two calculations along the lines 
described above as effective medium approach will lead us to the solution of the full 
problem posed by Hamiltonian (0). 



Influence of spin-flip scattering 



in a two-band Hubbard model 



6 



2.2. Spectral density approach 

According to the effective medium ansatz, the self-energy part S^^^ ^(i?) connected with 
the intraband Hubbard interaction will be calculated using the Hamiltonian ([10|) . We use 
the spectral density approach (SDA) to solve this problem. This method is numerically 
simple and fast, and, as shown in the limit of infinite dimensions, its magnetic properties 



resemble qualitatively the quasi-exact Quantum Monte Carlo calculation |]T8| ^ 
The SDA has been studied for the single-band Hubbard model extensively [^, |3^ . 
In the following we will give only a short outline of the calculations. Starting point is 
the single-electron spectral density being defined by 

+00 

" —00 
where [...,...]+ denotes the anticommutator and (. . .) the thermodynamic average. The 
construction operators are taken to be in the Heisenberg time-dependent picture. 

In an exact spectral-moment analysis in the limit U ^ 00, Harris and Lange have 
shown that the spectral density essentially consists of a two-peak structure [^. Since 



ferromagnetism is widely believed to be a strong coupling phenomenon, any reasonable 
approximation aiming at ferromagnetism should contain this limiting case [Q. 

In the SDA, one makes the following ansatz for the spectral density, which will turn 
out to correctly reproduce the positions and weights of the quasiparticle peaks according 
to the Harris- and Lange-calculation in the limit of ?7 -^00. 

Sk,raAE) = Yl - E^^mJ (16) 

The unknown parameters E^^]^ ^ and o-; the quasiparticle energy and spectral weight 
can be calculated by the moment method. That means they are fitted by the use of the 
first four moments of the spectral density, which represent several sum rules, and which 
can be calculated directly from the Hamiltonian. 

+00 

m£,,,= j dEE-S^^^AE) (IT) 
—00 

= ([[... K„,.,i/]_,...,//]_,a^J+) (18) 
V ' 

n-fold commutator 

This procedure is identical to the one performed in ||26| for the conventional Hubbard 



problem. An explicit description of the calculation is presented there. As a result one 
obtains a self-energy of the following structure: 

The decisive terms are B^-a and Fk^m~a which distinguish this self-energy from the 
Hubbard-I solution p. There these terms would be replaced simply by the center 
of gravity of the appropriate Bloch band. Bm~a and Fk^m,-a niainly consist of higher 
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correlation functions. They may provoke a spin-dependent shift and/or deformation of 
the bands and may therefore be responsible for the existence of spontaneous magnetism 



34, IC]. The fc-dependent term -Ffc^.-cr seems to be of minor importance for the 



magnetic behaviour |]26|. Since X]fc-^fc,m,-o- = it does not change the center of gravity 
of the density of states. It is mainly responsible for a deformation and narrowing of the 
bands. We have therefore neglected this term in the following calculations. The term 
Bm-a has the following structure: 

Bm,a = J^^i^ exp{-ik{Ri - Rj))era,a{.k) J (4m,f7«j,m,a(2 Ui^rn-a " 1)) (20) 
i,j \ k / 

with em,o-(fc) = tm{k) + '^^^In ai'^i) ■ Fortuuately, this two-particle correlation function is 
accessible via the single-electron spectral density, no higher Green's functions have to 



be calculated |41]. One obtains the following expression: 



Bm,a ^ 



^^jW(^ / 2 X _ \ 

— oo 

This leads to a set of equations which can be solved self-consistently. Despite its 
obvious restrictions, e.g. the complete neglect of quasiparticle damping, the two-pole 
approximation together with the moment method is able to describe the magnetic 
properties of the Hubbard model surprisingly well |3^. Since the subject of this 
paper is the influence of the spin-flip processes on the ferromagnetism introduced by the 
on-site intraband Hubbard interaction, this choice of a numerically simple procedure 
here seems reasonable. However, one should bear in mind that conceptually, the SDA is 
a strong coupling method which surely becomes questionable for intermediate to weak 
couplings. Consequently, we restrict all the following considerations to situations with 
Um substantially larger than the free bandwidth Wm- 



2.3. Rigid Spin Approximation 

Next, we have to solve the Hamiltonian (pi]). The approximation scheme for this 
effective problem has to be chosen very carefully since our investigation aims at effects 
directly induced by the interaction (^). There is no standard method to solve this model 
beyond mean-field level. In the following, we want to apply a non-perturbative, moment- 
conserving, self-consistent method which explicitly includes spin exchange scattering. 

The basis of the approximation scheme is the similarity between the interaction (^) 
and the well-known sf- or Kondo-lattice model. The difference of the two models lies 
in the electron spin operator cTj In our model, this operator is built from electron 
construction operators. In the Kondo model, the charge degrees of freedom of the /- 
spin have been projected out. Only a pure spin operator remains. A formal equivalency 
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between both models can be reached if one artificially fixes the operator (Ti fn to its spin 
degrees of freedom. This can be done in an iterative way for both bands separately. 
That means for calculating the self-energy for band m one has to fix the spin-operator 
of the other band (cTj fjj) and vica versa. Thus no constraints are introduced concerning 
the quantum mechanical attributes of the partial electron system under consideration, 
as e.g. the non-distinguishability of the particles. In our opinion, the name Rigid Spin 
Approximation (RSA) is an appropriate denotation for this method. 

Now we can apply a proper approximation developed for the s/-model. We 
use a moment-conserving, self-consistent interpolating equation-of-motion decoupling 
scheme |]35| , ^ A discussion of the neccessary approximations and their implications 



can be found there, so that we can restrict ourselves to a short summary in this 
paper. In the equation of motion for the single-electron Green's function Gij^rn,a{E) = 



I \ k 



— 2J{Ei,i,j,m,u{E) + z^Ti^ij^jn^E)) (22) 

two higher Green's functions are introduced, the "spin-fiip" function Fi^ij^rn,a{E) = 
'i^i,m(^i,m-a;a]^^A) and the "Ising" function Ti^i^j^rn,<j{E) = {{<^lm(^i,my,aj,m,a))- 
these two Green's functions, the respective equations of motion can be obtained without 
problems. In each of them, new Green's functions are introduced. For these we 
carefully apply a sophisticated decoupling scheme. The central idea is to express all 
local higher Green's function in the equations of motion of the spin-fiip and Ising 
function by interpolation "ansatzes" between several non-trivial limiting cases, such as 
5" = |, ferromagnetic saturation of one band, empty and full band. The corresponding 
interpolation parameters can be obtained using a moment method similar to the one 
described in section |2]^. Within this approximation scheme we obtain a selfenergy of 
the following structure: 

^iiliE) = J{al^) + j'H^^^liE), (n„,.), {a+^arj, . . .) (23) 

The first term corresponds to the mean-field solution. The second part corresponds 
to higher order terms in J. The complex functional depends on several correlation 
functions, as e.g. the interband spin-exchange correlation (o"j^mCr~^) or the interband 
Ising correlation (cfj^af^), the selfenergy part Sm.^-E") and of course, via the effective 
medium, on the intraband selfenergy part T!"m}a{E). All correlation functions as well as 
the selfenergy parts have to be determined selfconsistently. The selfenergy part Sm,L(-E) 
is fc-independent due to the neglect of magnonic excitation energies, which are small 
compared to the electronic excitations under consideration [^. It is worth to mention 



that this method can be continued smoothly to the exactly solvable non-trivial limiting 
case of one electron in a ferromagnetic saturated background of /-spins. 
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3. Results and Discussion 



In this section we present the results obtained by the theory described above. The results 
are compared with the usual mean-field calculation for the spin-exchange interaction. 
The two-band model under consideration consists of two bands m = {0,1} of unit 
width (VF{o,i} = 1.0) thus defining the energy unit used in this paper. The bands are 
not degenerate, the centers of gravity are shifted by 0.1. As an example we choose 
tight-binding bcc free densities of states [^. In the single-band Hubbard model, the 



existence of ferromagnetism depends on the lattice structure, as has been shown by 
various methods (stability of the Nagaoka state |l^, SDA QMC in = oo |^3|). 
From these investigations, it follows that the system has stronger tendency towards 
ferromagnetism in non-bipartite lattices (e.g. fee). It seems reasonable to believe that 
this will also hold for the two-band model although we did no systematic analysis of the 
lattice dependence. 
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Figure 1. Magnetization m as function of band occupation n of a single-band Hubbard 
model with [/ = 5, T = on a bcc-lattice. The dashed line indicates saturation 
(to = n). Antiferromagnetic ordering is expected close to half filling {n ~ 1) which 
is not considered in this paper. Due to particlc-holc-symmctry, it is sufficient to 
investigate the region of < rt < 1 



One expects the local intraband Coulomb matrix element U to be large whereas the 
interband exchange coupling constant J could be one order of magnitude smaller. As 
mentioned above the SDA which we use to find the intraband selfenergy part sjf2i^(£') 
is basically a strong-coupling theory. In the large-U regime, however, magnetic key 
quantities such as the Curie-temperature, the T = 0-moment, and so on, are already 
saturated, i.e. no longer ?7-dependent (see e.g. figure 11 in ||2^). More interesting 
is the J-dependence. We therefore restrict ourselves to the representative value of 
U = 5 and inspect in detail the infiuence of the interband interaction J. Furthermore, 
the evaluation is confined to T = 0, although the theory of course holds for finite 
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temperatures, too. 



3.1. The J = case 

In the case of vanishing interband couphng, J = 0, the situation is identical to two 
separate single-band Hubbard models which are only coupled by a common Fermi 
energy. It determines the respective partial band occupations ?t.{o,i} according to the 
total number of electrons per site ntot = uq -\- ui. For the single-band Hubbard model, 
the existence of ferromagnetism has been since long a matter of controversial discussions. 



but recent results confirm its stability for certain parameter regimes [|T5|, 0. The SDA 
has turned out to be able to reproduce the QMC results in this limiting case on a 



qualitative level [34, 38 



First, let us remind of a result obtained within a single-band Hubbard model. 
In figure |T], the dependence of the magnetization on the electron density (occupation 
number) is plotted. For n < n^c^^^^ ^ 0.56 the system is paramagnetic, only for 
n > n^c'^^^ ferromagnetic ordering is possible. With increasing n, the system becomes 
quickly saturated. Antiferromagnetic ordering occurs only in the very vicinity of half 
filling (n = 1.0) which we will never consider in the following discussion. More 
information on ferromagnetism in the single-band Hubbard model can be found in 

ii^ 



i9|, ^ ^ ^, y , |4| . 
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Figure 2. Magnetization m as function of total occupation ntot of a two-band system 
(solid line: total magnetization, dashed line: polarization of the lower band (m = 0), 
chained line: polarization of the upper band (to = 1)). The thin lines represent the 
partial occupation numbers per band. System parameters: U — b^J — Q,T — Q 



In figure 0, the magnetization as function of the electron density is plotted for the 
two-band situation. Together with the total magnetization mtot the partial occupation 
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numbers and polarizations per band are shown. The interband interaction is still set 
to zero (J = 0). As already mentioned above, both bands are coupled only via the 
chemical potential. Ferromagnetic order sets in in band m = 0, the lower band, when 
this band reaches a critical occupation of uq = n^f^^^^ ~ 0.63 > n^c^^^^ which is the 
case for a total occupation of utot ~ 0.98. The critical electron densisty of the lower 
band n^^^'^^ is larger than in the single-band case (n^^^^^). Even though the upper 
band is neither explicitely coupled to the lower band nor is ferromagnetic by itself, 
it has an influence on the para- to ferromagnetic transition of the lower band. This 
can be understood, since the upper band can act as source or sink of electrons for 
the lower band. With the onset of ferromagnetism, a spin-dependent band splitting 
takes place. This provides for a rearrangement of electrons between the bands due 
to the same chemical potential for both bands. Now, if the system tries to order 
ferromagnetically, when band m = reaches the "single-band critical occupation" of 
approximately Uq ~ 0.56, the corresponding shift of the densities of states will make the 
ferromagnetic phase instable. This happens until the real critical value of the two-band 



situation, uq = ^ 0.63 is reached. This effect is analogous to the findings in |45 

where a similar situation was examined using a Stoner-like theory. 

3.2. The Utot -dependency of the magnetization 



Meanfield-Approximation Rigid-Spin-Approximation 
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Figure 3. The same as figure ^ but witli finite J = 0.2. Left-hand side: Interband 
exchange interaction in mean-field approximation, right-hand side in RSA 



Let us focus on the interband-exchange interaction. Figure |^ shows, similar to 
figure H, the magnetization as function of electron occupation but with finite J = 0.2. 
On the left hand side, calculations were made using the mean-field approximation for the 
interband-exchange self-energy part, on the right hand side the rigid spin approximation 
as explained in chapter |2]1 was applied. 
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The mean-field result is easily understood: When the lower band reaches a critical 
occupation number uq = n^^^ ~ 0.50 < n[^(^'^^ a transition to a ferromagnetic state 
occurs. Now the interband exchange in mean-field approximation will provide for a rigid 
spin-dependent shift of the m = 1 quasiparticle bands proportional to the magnetization 
of the lower band. This induces a polarization of the upper band as well. The same 
mechanism now works as a feedback onto the lower band, increasing the magnetization 
even more. The lower band becomes thus quickly saturated. 

When analyzing the same situation using the rigid spin approximation, the picture 
gets modified in a rather drastic way. The onset of ferromagnetism is indicated by a 
critical occupation for the lower band n^c^^^^ which is only a little larger than nc™^"*. But 
the polarization of the upper band is much weaker than in the mean-field case. And 
furthermore, the lower band never reaches saturation. The latter is not a result of the 
weak polarization of the m = 1 band, but both effects have the same origin. Continuous 
spin-flip scattering prohibits the lower band from reaching saturation. These processes, 
in addition to the generally stronger quasiparticle damping also reduce the magnetic 
polarization of the upper band. 



3.3. The J -dependence of the magnetization 



0.70 




Figure 4. Total magnetization m (solid line) and partial polarization of the two 
bands (dashed: lower band (0), chained: upper band (1)) as function of the interband 
coupling constant J for fixed total occupation ntot = 1-0 and U = 5, T = 0. 



Next, we want to investigate the J-dependence of the magnetization. This turns out 
to yield unexpected, non-trivial results. Figure ^ shows the magnetization as function 
of J for fixed [/ = 5 and utot = 1.0. In figure ^ the behaviour for very small J is 
plotted together with results obtained using mean-field approximation for the interband- 
exchange, and finally, the corresponding quasiparticle densities of states are plotted in 
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figure ^ For the chosen total occupation number = 1.0, the lower band is already 
ferro magnetically ordered for J = 0, whereas the upper band is still paramagnetic (see 
figure I). 

The magnetization curve in figure § can be separated roughly into 3 regions. For 
very small J the magnetization rises (region A), for intermediate J it decreases (region 
B) and shows a reentrant behaviour for the largest J under consideration (region C). 

3.3.1. Region A: The first of these regions is characterized by an increasing 
magnetization as seen more clearly in figure |^. There, additionally to the RSA result, 
the mean-field curves are plotted. These are simply understood, any finite J will induce 
magnetic polarization into the upper band, and via a feedback, push the lower band 
into saturation. Band m = 1 will become more and more polarized, finally reaching 
saturation, too (not plotted in figure |^). The same effect can be seen in the RSA 
calculation, which is in first order of J identical to the mean- field result. 

3.3.2. Region B: Already for J ^ 0.05, which is still a very small parameter, deviations 
from the mean-field result are quite strong. This leads to the second regime in figure 
^. Here, the magnetization decreases, first slightly, but with increasing J stronger and 
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Figure 6. Quasiparticle densities of states (QDOS) for certain values of J for the 
same system parameters as in figure 3 with the chemical potenial /i = 0. The m = 
bands are drawn in solid, the m = 1 quasiparticle bands in dashed lines. Spin 1 and 
spin I QDOS are plotted separately as indicated. 
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stronger. In this parameter regime, the higher-order contributions in J to the self- 
energy become more and more important. The connection between the reduction of 
the magnetization and spin-flip processes becomes obvious by a comparison with the 
conventional s/-model. The important observation is the onset of a small dip in the 
quasiparticle densities of states (QDOS) with increasing J (figure |). This dip is on a 
much smaller energy scale than the Hubbard splitting, which is also clearly visible in 
figure |. 

To understand this feature, we want to relate it to a feature known from the 
conventional s/-model. There, a bandsplitting can occur for intermediate to large values 
of J. The size of the gap scales roughly with J. The physics causing this gap can be 
understood best by examining the exactly solvable special case of one electron in a 



saturated spin background |4g, |36|. In this case two different elementary excitations 
can be observed. One represents the scattering of an electron accompanied by the 
emission of a magnon, whereas the other can be connected with a bound state of an 
electron with a cloud of magnons. The latter manifests itself in the spectral density 
by a delta-like peak splitting of the scattering part for large enough values of J. The 



respective quasiparticle is called magnetic polaron |3^. For the general case of the sf 



model, i.e. for finite electron density and a not fully polarized spin system, a similar 
band splitting due to the same two elementary processes will occur for large enough 
values of J in both spin-resolved sub-bands |^. The dip seen in the QDOS of figure ^ 
is a precursor to the magnetic-polaron-induced bandsplitting. 

The reduction of the magnetization originates clearly from the spin-flip terms in 
the Hamiltonian. A further test is the artificial neglect of quasiparticle damping by 
setting Q'Ecr(-E') = 0. Even then, a reduction, though smaller, of the magnetization with 
increasing J is found. So additionally to the generic damping effects, explicit spin-flip 
scattering depolarizes the system. 

3.3.3. Region C: Finally, one observes a reentrant behaviour of the ferromagnetic 
ordering in figure ^, the origin of which is not completely clear. One possible mechanism 
for supporting magnetic order would be a RKKY-like interaction mediated by the 
interband interaction (^. The reentrant behaviour would indicate, that only for 
J > 0.45, the effective RKKY-interaction is strong enough to have an effect. To support 
our proposal, we have performed a calculation for the case of = and finite J. Even 
for these more or less unphysical parameters, we find stable ferromagnetic solutions. 
But as in the finite-[/ case, these only occur for relatively large values of J > 0.64. So 
even this artificially castrated model without direkt Coulomb interaction shows similar 
behaviour, which is, of course, quantitatively modified in the more realistic model with 
finite U . 

How this proposed mechanism is related to an enhanced stability of ferromagnetism due 



to the two-band situation found in ||2^ is an interesting, but open question. 
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4. Conclusions 

In this paper, we have analyzed a special multi-band Hubbard model which contains 
the two interaction terms believed to have the biggest influence on the stability of 
ferromagnetic ordering. 

Our model (|l]) is clearly insufficient to describe the rich physics of transition 
metals. But the same applies probably to any other model Hamiltonian which is still 
tractable within many-body theory. However, ab-initio calculations as e.g. density 
functional theory applied within the local density approximation (LDA) appear to 
underestimate just those correlation effects which seem to be decisive for phenomena 
like ferromagnetism. A proper combination of LDA calculations with a many-body 
treatment gives a promising way to solve these difficulties. The LDA calculation 
accounts for all interactions on a mean-field level, the many-body treatment should 
restrict itself only to the most important correlations. We believe that the model (|l|) 
contains in this sense those interactions whose contributions beyond meanfield have the 
biggest impact on magnetism. A combination of LDA calculations with the model (|l|) 
using a simpler approximation on the exchange part yields very good agreement with 
experimental facts 0, |^, ^ 0. For example, by fitting the interaction constants U and 
J to groundstate properties, the Curie-temperature could be calculated astonishingly 
accurate. 

In this paper, we introduced a more sophisticated approximation scheme, which is 
especially improved on the exchange interaction part of model (|I]). Our analysis was 
based on the fact that the intraband Hubbard interaction alone is able to form ordered 
magnetic moments in a band (itinerant magnetism). We investigated the infiuence of 
interband exchange coupling, often referred to as Hund's rule coupling, on the stability 
of ferromagnetism. Only for very small values of J, the respective interband coupling 
constant, we can verify the mean- field result which leads always to an enhanced stability 
of spontaneous ferromagnetism. Already for J ^ 0.1 which we call intermediate coupling 
the magnetization gets suppressed by spin-fiip scattering (using the free bandwidth as 
energy scale: W = 1). With increasing J the ferromagnetic order almost vanishes. The 
reentrant behaviour found for even larger J might be due to an RKKY-like ordering 
mechanism mediated by the interband interaction (^. The regime of J most often 
referred to transition metals belongs to intermediate coupling strengths. Our results 
indicate that the infiuence of the spin fiip processes can manifest itself in a rather 
dramatic reduction of the magnetization. 
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